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GIBBS MEASURES AND FREE ENERGIES OF ISING-VANNIMENUS 
MODEL ON THE CAYLEY TREE 

FARRUKH MUKHAMEDOV AND HASAN AKIN 


Abstract. In this paper, we consider Ising-Vannimenus model on a Cayley tree for order two 
with competing nearest-neighbor, prolonged next-nearest neighbor interactions. We stress that the 
mentioned model was investigated only numerically, without rigorous (mathematical) proofs. One 
of the main point of this paper is to propose a measure-theoretical approach the considered model. 
We find certain conditions for the existence of Gibbs measures corresponding to the model. Then 
we establish the existence of the phase transition. Moreover, the free energies of the found Gibbs 
measures are calculated. 
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1. INTRODUCTION 

Gibbs measure is ne of the central objects of equilibrium statistical mechanic, a branch of prob¬ 
ability theory that takes its origin from Boltzmann [2], Also, one of the main problems of the 
Statistical Physics is to describe all Gibbs measures corresponding to the given Hamiltonian. It is 
well known that such measures form a nonempty convex compact subset in the set of all proba¬ 
bilistic measures. The purpose of this paper is to investigate Gibbs measures of the Ising model 
[39j with ternary prolonged and nearest neighbor interactions on Bethe lattice of order two and to 
describe its extreme elements (pure phases). In [16] . we have studied phase diagram and extreme 
Gibbs measures of the Ising model on a Cayley tree in the presence of competing binary and ternary 
interactions. In [IS], we have obtained the extreme Gibbs measures of the Ising-Vannimenus (IV) 
model [39] without proof a number of theorems related to the description of the general structure 
of extreme Gibbs distributions of the Ising model on a Cayley, we studied Gibbs states (phases) 
that correspond in probability theory to what are called Markov chains with memory length 2 by 
using the method in [18]. In this paper, we combine the results obtained in [TT] and [18]. In [32] . 
the authors study the existence of p-adic quasi Gibbs measures by means of investigation of the 
obtained equations via methods of p-adic analysis. Note that the methods used in [32] are not valid 
in a real setting. Therefore, in this paper we will give a rigorous description of Gibbs measures 
corresponding to Ising-Vannimnus in the real setting. 

Recently, Ganfoldo et al [21] have obtained some explicit formulae of the free energies (and 
entropies) according to boundary conditions (b.c.) for the Ising model on the Cayley tree. Also, 
Ganfoldo et al m study the Ising model on a Cayley tree. They show that a wide class of new 
extreme Gibbs states is exhibited. 

Until now, many researchers have investigated Gibbs measure with memory one over Cayley 
tree. In this paper, we consider Ising model on a Cayley tree for order two with competing nearest- 
neighbor, prolonged next-nearest neighbor interactions. We stress that the mentioned model was 
investigated only numerically, without rigorous (mathematical) proofs. We propose a rigorous 
measure-theoretical approach to investigate a Ising-Vanniminus model on a Cayley tree of order 
two. We find certain conditions for the existence of Gibbs measures corresponding to the model. 
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Then we establish the existence of the phase transition. We investigate the behavior of free energy 
at these new transition points. We present some explicit formulae of the free energies corresponding 
to boundary conditions for the Ising-Vannimenus model on the Cayley tree of order two. 


2. Preliminaries 

2.1. Vannimenus Model with competing interactions on Cayley tree. A Cayley tree T k 
of order k > 1 is an infinite tree, i.e., a graph without cycles with exactly k + 1 edges issuing 
from each vertex. Let us denote the Cayley tree as T k = (V, A, i) where V is the set of vertices 
of T k , A is the set of edges, and i the incidence function which associates with each edge l £ L 
its endpoints. Two vertices x and y,x,y£V are called nearest-neighbors if there exists an edge 
l £ A connecting them, which is denoted by l =< x,y >. The distance d(x,y),x,y £ V, on the 
Cayley tree r fc , is the number of edges in the shortest path from x to y. Two vertices x,y £ V° 
are called the next-nearest-neighbours if d(x, y) = 2. Next-nearest-neighbour vertices x and y 
are called prolonged next-nearest-neighbours if they belong to the same branch which is denoted 
by > x,y <. Let spin variables cr(x),x £ V, take values ±1. Then the Vannimenus model with 
competing nearest-neighbours and next-nearest-neighbours binary interactions is defined by the 
following Hamiltonian 

(2.1) H ( a ) = — J p ^2 a ( x ) a ( y ) ~ J ^2 cr ( x ) c7 (y)’ 

>xiy< <x ’y> 

where the sum in the first term ranges all prolonged next-nearest-neighbours and the sum in the 
second term ranges all nearest-neighbours. Here J pi J £ R are coupling constants. 

The distance d(x, y),x, y £ V, on the Cayley tree T k , is the number of edges in the shortest path 
from x to y. For a fixed x° £ V we set 

W n = {x £ V\d(x, x°) = n}, V n = {x £ V\d(x, x°) < n} 

and L n denotes the set of edges in V n . The fixed vertex x° is called the 0-th level and the vertices 
in W n are called the n-th level and for x £ W n let 

S(x) = {y £ W n -\\d(x, y) = 1}, 

be the set of direct successors of x. 

Note in [39] it is assumed that J > 0 and J p < 0. Below we consider model (1) with arbitrary 
sign of the coupling constants. 

As usual, one can introduce the notions of Gibbs distribution of this model, limiting Gibbs 
distribution, pure phase ( extreme Gibbs distribution ), etc (see m, mm , ra- 


3. Translation-invariant Gibbs measures of the Ising-Vannimenus model 


In this section we define a notion of Gibbs measure corresponding to the Ising-Vannimenus model 
in a general setting, i.e. for arbitrary nearest-neighbor models (see [39]). We propose a new kind 
of construction of Gibbs measures corresponding to the Ising-Vannimenus model. 

Let <I> = {—1, +1} (3> is called a state space) and is assigned to the vertices of the tree = (V, A). 
A configuration er on V is then defined as a function x £ V —>• <x(x) £ <&; in a similar manner one 
defines configurations a n and ui on V n and W n , respectively. The set of all configurations on V 
(resp. V n , W n ) coincides with H (resp. Qv n = <& Vn , ^w n = d> l|/n ). One can see that 

Qy n = lV n -i x Qw n - Using this, for given configurations <7 n _i £ and u £ &w n we define 

their concatenations by 


(<7 n _i V w)(x) 


On-l(x), if X£V n -l, 
w(x), if x € W n . 
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It is clear that cr n _i Vw€ ^v n - 

For the Ising model with spin values in = {—1, +1}, the relevant Hamiltonian with competing 
nearest-neighbor, prolonged next-nearest-neighbor and two-level triple interactions has the form 

(3.1) H(a) = —J p Y a(x)a(y) — J Y a(x)cr(y), 

>x^f< <x,y> 

Assume that h : V \ M is a mapping, i.e. 

^■xy,a(x)a(y) = {hxy,++> h xyr \—, h X y t — \-,h X y t ), 

where h. xy ,a{ x )a( y ) € M (a(x),a(y) G {-1,+1}) and x,y G V\ {z (0) }. 

Now, we define Gibbs measure with memory of length 2 over Cayley tree. 

(3-2) = ^-exp[-/3H n (a) + Y Y cr(x)a(y)h xy ^ x) ^ y )]. 

n xGW n -i y£S(x) 

Here, as before, f3 = -^ and a n : x G V n —> cr n {x ) and Z n is the corresponding to partition 
function 

Z n= Y ex P [~/3H(cr n )+ Y Y a ( X ) a (.y) h xy,*(xMy)}- 

CTn&^Vn xeWn -1 y&S( X ) 

In this paper, we are interested in a construction of an infinite volume distribution with given 
finite-dimensional distributions. More exactly, we would like to find a probability measure // on H 

(n) 

which is compatible with given ones , i.e. 

(3.3) y,(a G H : a\ Vn = cr n ) = ™\ a n), for all a n G n G N. 

The consistency condition for y. n ,h{&n)i n > 1 is 

(3-4) Y ^K-iVwl^/f^K-i), 

ueQwn 

for any cr n _ 1 G ^v„-i- 

This condition according to the theorem implies the existence of a unique measure defined on 
H with a required condition (13.311 . Such a measure is said to be Gibbs measure corresponding 
to the model. Note that more general theory of measures has been developed in [23]. 

For cr n _i G V n -\ and rj G W n , we can define the Hamiltonian as following; 

(3.5) H n (a n -1 V rj) = -J Y a ( x ) a (y) ~ J Y Y a ( x ) r l(y) 

<x,y>eVn -1 xeWn-i y€S(x) 

- j p y cr (®)°'(i/) - j p Y Y a ( x ) r i( z ) 

>X,y<€Vn-l xGWn -2 zeS 2 (x) 

= H n (a n - 1) - J Y Y °'( x ) 7 ?(y) - J p Y Y 

xeWn-i yeS(x) xeWn -2 zeS 2 (x) 

Lemma 3.1. If ^ g = At and g = then there exists D G R such that a = DN\, 

b = DN 2 , c = DN 3 and d = DN 4 . 

Proof. Due to statement of Lemma, we have 
a = DNi, b = DN 2 , 
a = DiiVi, c = DiN 3 , 

a = D 2 Ni, d = D 2 N 4 . Therefore, from § = D = D x . Similarly, from § = = Al, 

D = D 2 . So, D = D\ = D 2 . □ 
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Theorem 3.2. Probability distributions n n (<j n ), n = 1,2,..., in (13.21) are compatible iff for any 
x, y £ V the following equations hold: 


(3.6) 


g hx y , + + + h x y , — + 


ghxy,++-\-hxy, + — 


n 

z€S(y) 


n 

zeS(y) 


exp[h y , ;++ ](gfr) 2 + exp[-h yZi+ _] 
exp[h y2i++ ]a 2 + exp[—h yz ,+-]b 2 

exp[—hy Z| —j_] + exp[h yz ,—](afr) 2 

exp[-h yz ,-+]b 2 + exp[h yZi —]a 2 


n 

z&S(y) 


exp[h yZi + + ](a6) 2 + exp[hy Z|+ _]] 
exp[—h yz -+]b 2 + exp[h :y2 ,—]a 2 ’ 


where a = exp (j3J) and b = exp (/3J P ). 


Proof. Necessity. From the equation (13.41) . we have 

L n exp [-/3H n (<r n _! V rj) + 0'(x')(j(y')}lxy,cr(x)cr(y) ] 

r]efiw n xeW n -1 yeS(x) 

= exp[—/3H n (a n -i) + ^ ^2 ^(x)cr{y)Ky,a(x)a(y)}, 

X£W n -2 yeS(x) 

where L n = E -1 . 
n z n 

From the equation (13.51) . we have 

(3.7) L n ^2 ex P[—pH n (a n -i) — fdJ ^2 ^ cr(x)r](y) 

y£t} Wn xGW n -i y£S(x) 

~ PJp E E cr{x)r]{z) + E E (j (x) c (y) h^y (T ( 3 ,) (T (y) ] 

xeWn-2 zgS 2 (x) xeW n -iyeS(x) 

= exp[-/3F n (cr n _i) + ^ (T ( x ) a (y) h xy,<T(x)<r(y)}, 

XGWn -2 y£S(x) 


n n n e exp[a(y)r](z)h yz ^ y)ri ^ + f5r](z)(Ja(y) + J p cr{x))\ 

xGWn—2 yeS(x) zeS(y) p(z)e{=Fl} 

= n n ex p{ (T ( x ) a (y) h xy,c7(xMy)}- 

x€W n -2 y&S{x) 


Let us fix < x,y > and considering all values of a(x),cr(y) € {—1,+1}. Then from (13.71) . we 
have 


(3.8) 


0 hxy,++-\-h'xy, — + _ 


E exp[? ? (^)(h y2i+ ^ (2) + /3(J + J p ))] 

y{z)e{zf 1}_ 

s , E exp[77(^)(h y2)+ ^ (2) + f3(J - J P ))} 

{y> p(z)e{=F 1 } 


n 


(3.9) 


0 hxy, - ~\~hxy ,-\— 


E ew[-y(z)(h yz _ v ( z) + P{J + J p ))} 

»7(z)S {=Fi}_ 

cQ , s E exp[-?7(^)(hy Z ,_„(*) - P(-J + 4))] 

zGb(y > p(2)G{T 1} 


n 


(3.10) 


E exp[ry(^)(h ?y2i+ ^ (2) + j3 (J + J p ))} 

y(z)e{- fl>_ 

^ q( , A E ex p[— ?y(^)(h y2 _^ (2) — (3(—J + J p ))} 

Z£b(y > p(2)€{ T l} 


n 
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These equations imply the desired equations. This completes the proof. 

Sufficiency: Now assume that the system of equations in (|3.6D is valid, then from Lemma 13.11 we 
get 


^ exp[a(y)ri(z)h yz ^ y) ^ z) +/3rj(z)(Ja(y) + J p <r{x))\. 

z£S(y) »?(z)6{=Fl} 

From the last equality, one can get 

( 3 .n) n n D{x : y s jg rT ^ x ^ cr ^^' x v 

xeW n - 2 yeS(i) 

e W{y)v{^yz,a(y) r i(z)+^{z){Ja{y)+J p a(x))}' 

xeWn -2 yeS(x) z€S(y) y(z)£{zf 1} 

Multiply both sides of the equation (13.111) by e ~P H n-i( a ) an d denoting 

Un= IJ II 

xGW„~ 2 y&S(x) 

we have from (13. lip . 



fiHn — 1 ( cr ) “E XI E CT ( ; r) cr (2/)^'xy,cr(x)f7(y) 

TT 1( o x61V„_2 yes(x) 


n n n 


l(«T) 


e 


,[o'(y)»7(z)h wz , CT (j ) )^( z )+/37j(z)(Jo-(i/)+Jp(r(a:))] 


xeVU„_2 y&S(x) z£S(y) 
Therefore, we get 


U n —lZn—lfJ"ii—l{ ( 7) — ^ ^ 6 
V 


V(z)£{ z El} 

-^JTi(erV?j)+ E E v{x)v(y)h xy ,<T(x)v(y) 

x£W n _2 s6S(i) 


(3.12) U n —\ Z n —\ y , n —\( y ( j^j — Z n YMu \/ y ) ■ 

v 

As // n (n > 1) is a probability measures, i.e. 

Mn-lM = ^ /X n (cr V 7?) = 1. 

cre{-i,+i } v "- 1 <re{-i,+i} v "-i ??e{-i,+i} w ™ 

From these equalities and the equation (13.121) we have Z n = U n -\Z n -\. 

This with the equation ()3.121) implies that (|3.4I) holds. □ 

According to Theorem 13.21 the problem of describing the Gibbs measures is reduced to the 
descriptions of the solutions of the functional equations (13.61) . 


4. The existence of phase transition 

Below we produce restrict ourselves to the case 

(4.1) h j X y^ —)—[— hxy t —y i and h X y : — h X y ,-\— ^ 2 ■ 

Denoting Since In u\ = h xyt++ = h xy , _f and In U 2 = h xy< _= h xy ^ _for any x, y € V, one can 

produce 


(4.2) 


u\ = 


a z b z u\U2 + 1 
a 2 uiU 2 + b 2 


= « 2 . 
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where a = e- and b = e~p. Hence u± = u z . Let u\ = U 2 = u, we have 


u = 


Let us consider a non-linear function: 

(4.3) f(u) = 


( ab) 2 ir + 1 
a 2 u 2 + b 2 

(■ ab) 2 u 2 + 1 
a 2 u 2 + 6 2 


Let us substitute as c = a 2 and d = b 2 , then we have a non-linear function as; 
(4.4) 


cdu 2 + 1 
9W = — 2 - 


cu 2 + d 

The analysis of the solution of the equations (13.611 is rather trick. In this subsection, we will 
study the translation-invariant solutions, where h xy = h is constant Vx, y € V. Let us consider 

cdu 2 + 1 

(4.5) u = — 5 --, 

v ’ cu 2 + d 

where u = e h . Note that if there is more than one positive solution for the equation (14.5p . then 
we have more than one translation-invariant Gibbs measure corresponding to the solution of the 
equation (14.51) . Fixed points u correspond to Markov chain Gibbs measures [41| (also splitting 
Gibbs measures or Bethe Gibbs measures). 

Proposition 4.1. The equation u = has one solution if either c < 1 or d < 3. If d > 3 

then there exist 771 (d), 172 (d) with 0 < 771 (d) < 772 (d) such that equation (14.51) has three solutions if 
771 (d) < c < 772 (d) and has two solutions if either 771 (d) = c or 772 (d) = c. 


Proof. Denote 


We have 


cdu 2 + 1 
9W = ~ZT2 - 


g\u) = 


cu 2 + d 
2 c (—1 + d)(l + d)u 


g"(u) = - 


(d + cu 2 ) 2 
2c(—1 + d)(l + d) (—d + 3cu 2 ) 


(d + cu 2 ) 3 

In particular, if d < 1, then g is decreasing and there can only be one solution of g(u) = u, 
therefore we can restrict ourselves to d > 1. We have that g is convex for 0 < u < and is 

concave for u > \J thus there are at most three solutions to g(u) = u. In fact one can see that 
there is more than one solution if and only if there is more than one solution to ug'(u ) = u, which 


is same as 


c 2 du 4 — c(d 2 — 3)u 2 + d = 0 

By means of the elementary analysis we get d 2 — 3 > 0 and A = c 2 (d 2 — 9)(d 2 — 1) > 0; and 


Ui = 


(d 2 - 3) + V(^ 2 -9)(d 2 - 1) 


2 cd 


u 0 = 


(d 2 — 3) — ^(d 2 — 9)(d 2 — 1) 


2 cd 


From the elementary analysis, one can see that then g'(u\) < 1 < g'(u, 2 ), if 


(-3 + 3d 2 + V9 - 10d 2 +d 4 J 


<c< 


(3 - 3d 2 + V9-10 d 2 + d*y 


32d 3 (d 2 - l ) 2 [d 2 - 3 + V9 - 10d 2 + d 4 ) ' 32d 3 (d 2 - l ) 2 [d 2 - 3- V9 - 10d 2 + d 4 ) 

In this case, we have three fixed points of the function g. □ 
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Figure 1. Jp = 34.6, J = -13., T = 27.5 



Figure 2. The region with three positive fixed points of the function (J4|) 

5. Free energy 

In this section, we study the free energy of depending on the boundary conditions for the Vanni- 
menus Ising model on Cayley tree. By previous sections we know that for any boundary condition 
satisfying the equations (13.6(1 there exist Gibbs measures corresponding to the IV-model. For this 
model the partition function is given by 

(5.1) Z n = Z n (/3,h)= ^2 exp {-/3H n (cr n ) + ^ ^ cr(x)cr{y)h xy ^ x) ^ y) }. 

crgn^ xGW n -1 y£S(x) 

Here the spin configurations a n belong to Q Vn and 

^ \h'xy,cr(x)a(y) € 2/ € F} 

is a collection of real numbers that stands for boundary conditions. In this section we will investigate 
the dependence with respect to boundary conditions of the free energy defined as the limit: 

F{P,h)= lim —j—In Z n (/3,h). 

rwoo fj\ V n \ 


(5.2) 
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Let us consider the equation 

D(x, y')e a ^ a ^ hxy ' !j( ' x '> a ( v '> = n E e uh yZ}Cr ^ y - )u +Pu(Ja(y)+Jp(7(x)) 

zeS{y) m€{±1} 

U n ~i = n n yi D ( x ’y)- 

x€W n -2 z£S(y) uG{±l} 

So, we have Z n = U n -\Z n -i. 

5.1. The free energies of translation-invariant Gibbs measures. In this section, we discuss 
the behavior of free energy as a function of the model in presence of an external field. We as before 
will consider the boundary conditions (14.11) . i.e. 

hxy,- 1 —h — ^xy, —T — ^l and h X y } — — hxy,-\ — — 1^2 ■ V <C x : y ^G L . 

Proposition 5.1. The free energies of compatible translation-invariant (TI) boundary condition 
exist and are given by 

(5.3) F TIl (/?, h) = — hi [cosh(/ij + j3(J + J p )) cosh(/i; + f3(J — J p ))], 

where hi are the varieties such that m = e hi are the fixed points corresponding to the equations 

m- 

Proof. From (14.11) one finds 

D(x,y)e hxy ’ ++ = n ^ff l yZt +++P(J+J p ) _j_ g— h yz ^ —— /3(J+J p )j 

z&S(y) 

(5.4) = I] cosh [ V ’++ ^ ^ + f)(J + J p) ). 

zeS(y) 


(5.5) 


D ( x , y)e hxy ,-+ = 11 j^g^y2,+++^(^ Jp ) _|_ g ^ yz , 2 — ^(4 4 p )j 

zeS(y) 

= JJ 2e hvz ,++ ~ 2 hvx ' + ~ cosh [ V++-V.+- 

zGS(y) 


+ /3(</ — «/p)] • 


Multiply the equations (15.41) and (|5.5D . then we have 
(5.6) 


D ( x , y ) = 4 b ( y , z ) 

zeS(y) 


where 


b ( y , z ) = (cosh [ '•»»,++ V+- +fK j + Jp) ] cosh [ '^++ V,+- + „(J _ Jp) ]) I. 


Un-l 


n n d ( x >v) 

xeWn-2 y£S(x) 

4 ,Wn -il JJ b ( x , y ) 

yeWn-1 z£S(y) 


4\w„-i\ 


E E lnb ( x , y ) 

3 V6V! 2 6 S(s) 
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where a(x, y ) = In b(x, y ). 

For this model, we can give the partition functions as follows; 


Z n — U n —\Z n —i 


(5.7) 


E E a(y,z) 


E 


— ^}Vn-l\ e yeW n _ 1 z£S(y) gSi6H'„_ 2 3 1 6S(s) 

d(x.y) 

— ^\Vn-l\ e <x,y>eV n 


E a(yi,zi) 


E E a(y,z) 

_ e y£W 0 z£S(y) 


Therefore, we get the free energy as follows; 

Frh (P, h) = - lim In D(x,y) 

n-^oo 

(5.8) = - In [e hl ~ h2 cosh(^-^-^ + /?(J + J p )) cosh(—+ /3(J - J p ))]. 

Assume that e hl = e hl . According to the formula (|5.8jl . for the case 1 of Vannimenus Ising model 
the free energies of translation-invariant (TI) boundary condition are given by 


(5.9) F TIl ((3,h) = ~ In [cosh(/ij + /3(J + J p ))cosh(^ + /3(J- J p ))], 

where hi are the varieties such that Ui = e hi are the fixed points corresponding to three cases. □ 


F 



Figure 3. J v = 34.6, J = -13, u\ = 0.260261, u 2 = 1.18483, u 3 = 3.52491 
The residual entropy at T = 0 is defined as follows: 


(5.10) 


Soo = ~ lim 

0—>co 


F{h) - F a 

j_ 

0 


where F \= lim F(h). 

0—>oo 

Let us compute the entropy 


(5.11) Stf,h) = - 


dF(/3, h) dF({3, h ) 1 


dT 


dp (3 2 


= — (J + J p ) tanh[/t + P(J + J p )] — (J — J p ) tanh [h + /3( J - J p )] 

The residual entropy at T = 0 is defined as follows: 
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(5.12) 


S 0 


— lim 

/3—>oo 


F(h ) - Foo 

i ’ 

13 


where Foo = lim F(h). 

/3—>-oo 

Let us compute the entropy 
dF(/3, h) 


(5.13)S(/3,Xi) = -- 


dT 


dF(0, h) 1 
dj3 0* 


In 9 In t • In T‘ 

= yr - (■ J + J v) tanh[^l + 0(J + J p )] + (J + J p ) tanh[^ - 0(J + J p )] 
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